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Abstract 

The random point field which describes the position distribution of the system of ideal 
boson gas in a state of Bose-Einstein condensation is obtained through the thermodynamic 
limit. The resulting point field is given by convolution of two independent point fields: 
the so called boson process whose generating functional is represented by inverse of the 
Fredholm determinant for an operator related to the heat operator and the point field 
whose generating functional is represented by a resolvent of the operator. The construction 
of the latter point field in an abstract formulation is also given. 



1 Introduction 

In the previous paper [TIJ, which we will refer as I, the authors gave a method which derives 
typical kinds of random point fields, the boson point process and the fermion point process 
on M. d , through the thermodynamic limit from random point fields of fixed finite numbers of 
points in bounded boxes in M. d . The purpose of the paper is to give the random point field which 
describes the position distribution of the system of ideal boson gas in a state of Bose-Einstein 
condensation [BEC] as an extension of I. 
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Let us consider the system of N free bosons in a box of finite volume V in IR d and the 
quantum statistical mechanical state for the system of a finite temperature. Regarding the 
square of the absolute value of the wave function as the distribution function of the positions 
of iV particles together with thermal average, we obtain a random point field of N points in 
the box. In I, the thermodynamic limit, N, V — > oo and N/V — > p, of the system for small 
p as well as the system of fermions for every positive p are taken to get the boson as well as 
the fermion point processes on Mr in a simple and straightforward way. As applications of the 
approach, the system of para-particles and the system of composite particles are studied. The 
argument is based on the unified formulation of boson/fermion processes of |STj . For general 
references of this field, see e.g. jBE] and references cited there in. 

In this paper, we study the case of large p (corresponding to BEC) which needs technically 
elaborate analysis for the largest eigenvalue go(L) of the deformed heat operator Gl in the box 
of size L and the saddle points zq and zq for complex integrals related to generalized Vere- 
Jones' formula |V| 1ST] more than those in I. As the result of the thermodynamic limit, we get 
a random point field on which are given by convolution of two independent point fields: 
1. the boson process whose generating functional is represented by inverse of the Fredholm 
determinant for an operator related to the heat operator; 2. the point field whose generating 
functional is represented by a resolvent of the operator. 

The paper organized as follows: In §2 the construction of the point field which appears in 
the resulting point fields as the second independent component (see above). The construction is 
made in a general framework of random point fields similar to jSTj . i.e., on the locally compact 
space of second countability. §3 devoted to the analysis of the thermodynamic limit in R d . 



2 Abstract formulation of the random point field 

Let R be a locally compact Hausdorff space with countable basis and A a positive Radon 
measure on R. We regard A as a measure on the Borel a-algebra B(R) which assigns finite 
values for compact sets. Relatively compact subsets of R will be called bounded. On L 2 (R; A), 
we consider a (possibly unbounded) non-negative self-adjoint operator K which satisfies: 

Condition K 

(i) [locally houndedness] For any bounded measurable subset A of R, the operator K x I 2 xa is 
bounded, where xa denotes the operator multiplying the indicator function xa- 

(ii) G = K(l + K)~ l has a non-negative integral kernel G(x,y) which satisfies 

G(x,y)X(dy)^l X-a.e.xeR. (2.1) 

For a measurable function / : R — > [0, oo) with compact support and a bounded measurable 
set A satisfying A D supp/, we have K 1 ^ 2 a/1 — e~? = A" 1 / 2 %av / 1 — e~f and hence that 

K f = {K l ' 2 ^i~e~fyK l l 2 ^l^e~l (2.2) 

is a bounded non-negative self-adjoint operator. Here we regard y/1 — e~$ the multiplication 
operator of the function expressed by the same symbol. Q(R) denotes the Polish space of all 
the locally finite non-negative integer valued Borel measures on R. 
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Theorem 2.1 For R,X and K which satisfy the above conditions and p > 0, there exists a 
unique Borel probability measure /ik, p on Q(R) such that 

[ e-<^dn K M) = exp ( - pWl-e-f, [1 + K f }~Wl ~ e~f)) (2.3) 

JQ{R) 

holds for any non-negative measurable function f on R with compact support, where ( • , • ) 
denotes the inner product of L 2 (R; A). 

Let us begin with some remarks before proving the theorem. It follows that G is self-adjoint 
and ^ G ^ 1, where 1 denotes the identity operator on L 2 (R; A). Without loss of generality, 
we may assume that the £>(i? 2 )-measurable function G(x,y) satisfies 

Vx,yeR: G(x,y)^0, G(x,y) = G(y,x) 

and 

VxeR: / G(x,y)X(dy) ^ 1. 
Jr 

Let us define the functions G n (x,y) inductively as 

G\x,y) =G(x,y), and G n+1 (x,y) = [ G n (x, z)G(z,y) X(dz) for neN. 

Jr 

Then we have 

Vx, y G R, Vn G N : G n (x, y) > 0, G n (x, y) = G n (y, x) 

and 

\/x G R, Vn G N : / G n (x, y) X(dy) < 1. 

It is obvious that G n (x,y) is the integral kernel of the operator G n for any n G N. 
Put 

n n 

K„ = ^G fe and K n {x,y) = Y,G\x,y). 

k=l k=l 

Then is the bounded non- negative self- adjoint operator which has non- negative integral 
kernel K n (x,y). The function 

oo 

K(x, y) = hm K n (x, y) = V y) (2.4) 

n~ >oo ^— ' 

fe=l 

is well defined, if we admit infinity as its value. 

Here we recall the following preliminary facts from functional analysis. 

Lemma 2.2 (i) Let 7i be a Hilbert space, C(J~L) the Banach space of all the bounded operators 
on H and {A n } ne fq a bounded increasing sequence of non-negative self-adjoint operators in 
C(7i). Then s-lim^oo A n exists and is a bounded non-negative self-adjoint operator. 

(ii) Suppose that A±, ■ • • , A n , • • • G C(L 2 (R; A)) converge to A G C(L 2 (R; A)) strongly, A n 
has integral kernel A n (x,y) for each n and 

< A n (x, y) t A(x, y) A® 2 - a.e. (x, y) G R 2 . 

Then A has A(x, y) as its integral kernel. 
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Proof : For (i), see e.g. |RJN| . 

For (ii), let / e L 2 (R;X). Then |/| e L 2 (i?;A) and (A n |/|)(x) = / A n (x, y)\f(y)\X(dy) 
holds. Taking the limit n — > oo (through a subsequence if necessary), we have (A|/|)(x) = 
/ A(x,y)\f(y)\\(dy) A — a.e. by strong convergence of the operators and the monotone con- 
vergence theorem. The a.e. boundedness of the integral in the righthand side ensures the 
identity for / instead of |/| by dominated (instead of monotone) convergence theorem. □ 

Now we have the following proposition. Here and hereafter, || • || and || • ||t stand for the 
operator norm and the trace norm for operators, respectively, and 1 1 • | \ p for the L^-norm for 
functions. 

Proposition 2.3 (i) Put K A = (K 1 / 2 x A )*K 1 / 2 x A f or bounded measurable A C R. Then, K A 
is a bounded non-negative self-adjoint operator and has K A (x,y) = XA(x)K(x,y)xA{y) as Us 
integral kernel. 

oo 

K A = Y,XAG k X A (2.5) 
fc=i 

holds in the sense of strong convergence of operators. 

(ii) For each k e N, H k = xaG((1 — XA)G) k ~ 1 XA is a bounded non-negative self-adjoint 
operator having non-negative kernel, denoted by Hk(x,y). R A = J2T=i^k exists in the strong 
convergence sense and is the bounded non-negative self-adjoint operator having non-negative 
kernel R A (x,y) = J2T=i H k( x >v)- 

(iii) R A = K A (l + K A )-\ \\R A \\ < 1. 

(iv) (1 + K A ) _1 Xa ^ a.e. holds, where we regard xa a s a function which belongs to 
L 2 (R;X). 

Remark : From (i) of the proposition and the argument above (|2.2|l . it follows that Kf = 
y/l — e~f K A y/l — e~f and its kernel is given by Vl — e~ J ^K(x } y)Vl — e - -^ for non-negative 
/ satisfying supp / C A. 

Proof : (i) Boundedness and self-adjointness of K A are obvious. 
Using the spectral decomposition K = J °° \dE\, we have 



Hence, 



yZxAG k XA\\= sup Y](XA0,G' fc XA0) 

sup / f t~~t) d (XA<P, E X XA<P) < sup f\ d(xA<t>, E x Xa<P) 

\K^x A ^\\\ = \\K A \\. 



sup 

1*112 = 1 



Since XAG k XA ^ holds for every k e N, Lemma l2~2T i) yields the existence of s-lim^oo J2k=i XAG k XA- 
On the other hand, thanks to the monotone convergence theorem, we get ()2.5|) in the weak sense: 



n f n / X \ k 

^XAG k XA(t>) = I (iTa) d( ^ x ^ ExX ^ 

k=l J k=l 
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Thus we have f|2.5J) in the strong sense. 

Lemma f2.2r ii) yields the assertion on the kernel of K A . 

(ii) It is obvious that H k is a bounded non-negative self-adjoint operator for every k G N. 
From the non- negativity of the kernel of G k , we have the non-negativity of the kernel H k (x,y) 
and 

< H k (x,y) ^ Xk{x)G\x,y) X K{y). 
Lemma f2.2lT j and the estimate 

n n „ 

Y^H k \\ = sup V/ ^jH k (x,y)<P(y)\® 2 (dxdy) 



k=l imia— ■ fc= i 



n „ n 

sup V/ |0(x)|xA(x)G fe (x )y )xA(y)|0(y)iA® 2 (dx^)^||y;xA^XAlKII^ 



we get the existence of the strong limit R A of {Y^=i Hk}n and its bounded self-adjointness. 
Lemma [2.2f ii) yields the assertion on the kernel of R\. 
(iii) From 

n n n 

Y,H k - J2xAG k XA = 5>aG[((1 - X A)G) k - 1 - G k ~ l ]xA 

k=l k=l k=l 

n k— 1 

= E E " XA)G) fe - z - 1 (-XAG)G'- 1 X A 

fe=2 i=l 

n—l n n—l n— I 

= ~ E E ^G((l - X A )G) fc - z - 1 X AXAG"xA = - E E ^-XaG'xa, 

Z=l fc=z+i i=i m=l 

we get the relation 

n n 

E f ) ~ E XA(x)G k (x, y)x\{y) 

k=l k=l 
n—l n—l 



EE / H rn(x,z)xh{z)G\z 1 y)xk{y)K dz ) a - e - 



1=1 m=l "'■ R 

in terms of kernels. Taking the limit n —>■ oo, we get 

R A (x,y) - K A (x,y) = - / #a(^, 2)^(2, y)K dz ) \ m - a.e.(x,y) 

JR 

by the monotone convergence theorem. It implies R A — K A = —R A K A and hence R A 
K A {1 + K^ 1 . Since K A is non-negative and bounded, ||-Ra|| < 1. 
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(iv) We may regard G as a contraction operator on L°°(R; A) because of (|2.1|) . H k is also 
contraction on L°°(R; A) for all fceN. Thus we have 



n n 

J2(H kX A)(x) ^ ^(tf feXA )(x) + { X aG((1 - X A)G) n -\l - xa))(x) 

k=l k=l 



n-l 

^ J2(H kX A)(x) + { X aG((1 - X A)G) n -\l - X a)) (x) < ■ ■ ■ ^ (xaG1)(x) ^ X a(x), 

k=l 

where non-negativity of the kernel of G and (j2.1)l have been used. On the other hand, we get 
Y^k=\ HkXA RaXa a -£- from (ii) through subsequence if necessary. Hence (1 + K a )~ 1 xa = 
Xa — RaXa ^ a.e. holds. 

□ 

(Proof of Theorem \2. ij) 

Recall that Kf = v 1 — e~$K A \J\ — erf, for non-negative measurable / and a bounded mea- 
surable set A D supp /. Since 

(1 + (1 - e~ f )K A )y/l ~ e~f(l + K f )~ Vl - e~f = 1 - e -/ = 1 - e~ f R A - e~ f (1 + Ka)' 1 
and 

1 + (1 - e~f)K A = (1 - e-fR A ){l + K A ), 

we get 



y/l-e~f(l + K f )~Wl-e-f = (1 + K A y\l - e^R^l - e~* R A - e~ f {l + K A )~ l ) 

oo 

= (l+KJ-^l-il-e-'R^e-fil+Kj,)- 1 ] = (l+K^-il+K^ 1 ^(^^^(I+^a)- 1 - 

n=0 

The Neumann expansion in the last step is valid since ||e~^i?A|| ^ II-^aII < 1- Hence we have 

-(y/l-e-f, [1 + K f ]~Wl-e-f) 

oo 

= -(xa, (1 + K a )- 1 Xa) + + K A )-\ A , e- f (R A e- f ) l (l + K A )- l XA ). 

1=0 

Substituting this identity to the right hand side of ()2.3jh expanding the exponential and sym- 
metrizing, we get a expression of the form 

°° 1 f 

V—/ a A n(x h --- ,x n )e~^^ s{Xk) dxi---dx n (2.6) 
t=t n] A* 

with a family of symmetric non-negative functions {0a«} for every A D supp /. For the existence 
of the measure [Ik, p on Q(R), it is enough to show the consistency condition [L]: 

°° 1 f 

a A n(xi, ■ ■ ■ ,x n ) = S2 - / cr {AuA)n+ i(x 1 , ■■■ ,x n ,y u ■■■,yi) dy x ■■ -dy h 
i=o 11 Ja > 



where A fl A = 0. This condition can be derived easily from the facts that the right hand side 
of ()2.3|) does not depend on A D supp / and that for a given A, {cta™} in ()2.6j) is uniquely 
determined a.e., since / can be arbitrary non-negative measurable function satisfying supp / C 
A. 

Thus we have proved Theorem 12.11 □ 



3 The Thermodynamic Limit 

In this section, we follow the arguments and the notations of I §2.2. However, let us review 
them briefly to make the article self-contained. 

Consider H L = L 2 (A L ) on A L = [-L/2,L/2] d C R d for d > 2 with the Lebesgue measure 
on A^. Let be the Laplacian under the periodic boundary condition in TIl- For k e Z d , 
<Pk ( x ) = L~ d l 2 exp(i27rA; • x/L) is an eigenfunction of Al, and { ip^' }k& d forms an complete 
orthonormal system [CONS] of TCl- The operator Gl = exp(/?Ai) has the integral kernel 



GUx,y) = ^e-*™My?\x)<pV{ y ), (3.1) 

for (3 > 0. We put gj^ = exp(— f3\2nk/L\ 2 ) which is the eigenvalue of Gl for the eigenfunction 
iff. L (x). We also need the operator G = exp(/3A) on L 2 (M. d ) and its integral kernel 

dp ^ m 2 +ip<x _ y) _ expj-jx -y\ 2 /4:f3) 



(27T) d (47T/?) rf / 2 

Let / : M. d — > [0, oo) be a continuous function of compact support. We will only consider the 
case where L is so large that Al contains supp/. We regard / as a function on Al naturally. 
Let 

G L = G L /2 e- f G L /2 , (3.2) 

where e~? represents the operator of multiplication by the function e~f . 

Suppose there are N identical particles which obey Bose-Einstein statistics in Al under 
the periodic boundary condition at inverse temperature (3. The basic postulates of quantum 
mechanics and of statistical mechanics of canonical ensembles yield 

Pl, n {xi, ■■■ ,x n ) = 7T^yP er {G(x u Xj)}^ j=1 (3.3) 

as the probability density distribution of the positions of N particles of the system, where Z B 
is the normalization constant and per represents the permanent of matrices. Here, we have 
set h 2 /2m = 1. We define the random point field ( the probability measure on Q(M. d ) ) //£ N 
induced by the map A^ 9 (xi, • • • , xjy) i— > J2f=i e Q(^ d ) from the probability measure on 
A^ which has the density 1)3. 3 J) . By Ef N , we denote the expectation with respect to /if N . The 
Laplace transform of the point process is given by 

J AN exp(- Y^f=i ffaflper {G L (xi, Xj)}ij=i dx 1 --- dxj 



J a n per {G L (xi, Xj)}^ j=1 dx x --- dx N 
f. N per {G L {xi, Xj)}f j=l dxi--- dx N 

(3.4) 



J AN per {G L (xi, Xj)}Fj =1 dx x --- dx 



N 
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Let us consider the thermodynamic limit, where N and the volume of the box tend to 
infinity in such a way that the densities tend to a positive finite value p: 

L, N -> oo, N/L d -> p > 0. (3.5) 

In this paper, we concentrate on the high density region 

f dp e-PW 2 , . 

where the Bose-Einstein condensation takes place. 

Theorem 3.1 (i) The operator K = G(l — G)^ 1 is a non-negative unbounded self-adjoint 
operator in L 2 (M. d ) and satisfies Condition K in §2. Moreover, Kf defined by \2.2\) is a trace 
class operator. 

(ii) The finite point fields defined above converge weakly to the random point field whose 
Laplace transform is given by 

w b r -</,£>! _ exp ( - (p - pJtyr^Ff, [1 + Kf^VT^T)) 

J " Det[l + K f ] ( ' 



in 



the thermodynamic limit \3. M'J.b]) . 



Remark: Thus the resulting point field of the theorem is a convolution of a point field which 
is an example of those discussed in §2 and a boson process. On the formulation of boson 
processes, we refer to |STj , where the operator K is assumed to be bounded, however the proof 
given there is also valid for the present case. 

Let us begin the proof with the following lemma, where we use the notation 



□l L) = -(k+(-K- 

k L\\ 2' 2 



for k G Z d . 



Lemma 3.2 For z G [0, 1], v — 1, 2 and L G [1, oo), let us define functions a u ( ■ ; z), a£ ( ■ ; z) 
on M. d by 

ze -P\p\ 2 

a ^ z ^= a- Z e-p\p\y 

and 



ai L \p;z) 

Then 



if P 6D[ L) 

a u (2irk/L; z) if p e D ( k L) for k e Z d - {0}. 



^ af\p- z) ^ a 1 {2p/{2 + Vd); 1) G L 1 (R d ) 
and the bounds for large L 

(c d (L/^P) d if d>4 

/ a { 2 L \p; z) dp ^ £(L) = I c 4 (L/VP) 4 log(cL/y% if d = 4 

U(VvW if d<A 



L" 



(27r) d 

hold, where Cd, c& and c are positive constants. 



8 



Proof: Since a u is monotone increasing in z and monotone decreasing as a function of \p\, we 
have 

a[ L \p-z) ^ sup ai L \p- 1) ^sup{a u ( q] l)\qem d ,L^l,\q\ > ^, \q - p\ ^ (27i/L)(Vd/2)} 



L>1 



^ sup{ a u (q; 1) | q G M d , L > 1, |g| ^ ^, |p| - ttv^/L < |g| }. 

In the case of \p\ ^ (2 + y/d)ir, the last supremum is attained at L = 1, \q\ = \p\ — n\fd then 
\q\ ^ |2p|/(2 + yd) holds. On the other hand, if \p\ < (2 + ^/d)ir, the supremum is attained at 
L = (2 + y/d)n/\p\, \q\ = 2n/L and then \q\ = \2p\/(2 + \fd) holds. For both cases, we get the 
bound ai\p- z) a v (2p/(2 + Vd); 1). Since d > 2, we get ai(2p/(2 + v^); 1) G L 1 (R d ). 
Integrating the angular variables, we have 

/ af\p-z)dp^ 1 ^ Vd [ a 2 (2p/(2 + Vd);l)dp 



(27T) 



I 7T -„2\r> dq — (- — 7= I Sdld, 



MVWJ J^/l (1 - e-«y \2txJW, 

where (3' = 4(3/(2 + Vd) 2 . Since 2 d ^ J^[q d - 1 e- q2 /(l - e~ q2 ) 2 }dq < oo for d > 4; l d < 
f^fljrplq*- 1 /^] dq = (4 - d)- 1 (L/ir^) i ~- d for d < 4 and 

J 4 ^ / 77 ~^Vi d 1+ f dq = const. + log— (3.8) 

we get the bounds for 7r-\//5 ^ L. □ 

(Proof of Theorem It is obvious that X = G(l — is a unbounded non- negative 

self-adjoint operator satisfying G = X(l + X) -1 . In fact, K is explicitly given by the Fourier 
transformation: 

K(f) = r^a^jl)^) 

for 

G DomfT = { ^ G L 2 (M d ) | a x ( • ; 1)^0 e £ 2 (M d ) }. 

Condition K(ii) for G are also obvious. 

Let us show the locally boundedness of K. For bounded measurable A C M. d , 

2 

oo 

2' 



\\^/K X k<P\\1 = \\^a~^l)HXK<P)\\l < llv^FTiyilall^Ca^ll 
^||a 1 (-;l)|| 1 || XA 0||^(27r) d p c || XA | 1211 - 112 
Thus -R' 1 / 2 xa is bounded. in (|2.4|) is given by 



oo oo „ 

K(x,y) = J2G n (x,y) = Y, / 

n=l n=l 



dp 



e -n^|p| 2 +ip-(x-j/) 
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(27r) d 

where we have used the dominated convergence theorem. From ai G L 1 (IR d ), K(x,y) is con- 
tinuous. The remark after Proposition 12.31 and the continuity of / yield that the kernel of 
Kf is continuous. Hence Kf is a trace class operator, because = J Kf(x,x) dx = 

p c \\l — e~^||i < oo. □ 

The rest of this section is devoted to the proof of the second part of the theorem. It is 
obvious from non-negativity of / that ^ Gl ^ Gl- Let us denote all the eigenvalues of Gl 
in decreasing order 

g Q {L)>~9i{L) >--->g,( L ) > ■ ■ 
Correspondingly, we relabel the eigenvalues { g£ }k& d of Gl as 

g (L) = l>g 1 (L)^---^g J (L)>--- . 

By the min-max principle, we have 

gj {L)>~ 9j {L) (j = 0,1,2,...). 

Note that has eigenvalue go(L) = g^ = 1. 
Put 

D L = G L -G L = Gf (1 - e-f)G)! 2 , W L = G x [ 2 ^\ - e~f, 



then D L = W L W* L . Note also that 

' / a£\p;z)dp = V Z9k = \ \zQoG L Qo(l-zQ G L Qo 



(2k) 



(3.9) 



Here Pq is the orthogonal projection on L 2 (Al) to its one dimensional subspace and 

Qo = l-Po- 
Now we have; 

Lemma 3.3 (i) \\Q G L Q (1 - Q G L Q )- 2 \\ T = £ fe ^ </*(£) /(l - g k {L)f < £(L) 

In the limitL — > oo, i/ie following convergences hold. 

(ii) L- d TrG L — > / Rd e-^l 2 dp/(27r) d = v 7 ^"", IPl||t — > ||1 - e"'| | 1 / v / 4^°' 
(hi) L- d ||g G L g (l - QoG L g )- 1 ||r = ^E^o^W/C 1 ~ flfcW) — Pe < oo 
(iv) // {z^} C (0, 1) and zl — > 1, £/ien 

su Px, y eA |[^l<5oG'l<5o(1 - ZxQoGxQo) -1 ]^,?/) - -> 

for any fixed bounded measurable set A C M. d . 



Proof : (i)-(iii) are immediate consequences of above remarks, Lemma 13.21 and the domi- 
nated convergence theorem. 
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For (iv), put e(p;x) = e ip ' x and 

e ( L \p;x) = e(2nk/L;x) if peOf' for k E Z d . 
Then Lemma f3. 21 and the dominated convergence theorem also yields 

/dp 
^p|e (L) (>; x-y)a[ L \p; z L )-e(p; x-y)ai(p; 1)| 

In the followings, we use the notation Bl = 0(L a ) which means 

3ci > c 2 > : ciL Q ^ 5 L > c 2 L a 

Lemma 3.4 

(i) For /ar#e L, g (£) - <jfo(£) 

= L~ d (Vl ~ e-f, [1 + W* L Q [1 - QoGLQol-'Qo^-Vl-e-/) (1 + o(l)) 
= (p< £) , (D £ - D L Q [1 - QoGlQoI^QoDl)^) (1 + o(l)). 

(ii) ||1 - e- / || 1 /i d = (^ L) , D^< L) ) > (^ L) ,I> L Qo[So(£) - QoGlQo^QoDup™) 



(iii) L d (o (L)-^ (L))G 



1 - e -/|| 1 (l + (l)) 
l + Pelll-e-Zllj 



II -e 



-/I 



(iv) Let (p Q be the normalized eigenfunction of Gl for eigenvalue go(L) such that 

,~(L) (L)x . n D , ~(L) (L) . / (L) /_(L) . _/ // (L) a n 

(^ L ) ? ^') = 0). r/iena = a' and \\<p'\\ 2 = \ \<p'\ \ 2 = 1 - a 2 = 0(L- 2 H(L)) hold. 



Proof: Here, we suppress the index L in gj(L),gj(L),(p and so on. First notice that 
(<Po,Dl<Po) = ||1 — e~f\\i/L d . From the min-max principle, d > 2 and the value of g± = 
exp(— P\2tt/L\ 2 ), we have 

9o = 1 > 9o > (<Po, G L ip ) = 1 - (<p , D L ip ) = l-6(L- d ) > 9l = 1-6(L- 2 ) > ~9i > ■ ■ • (3.10) 

for L large enough. Hence the eigenspace of Gl for its largest eigenvalue go is one-dimensional. 
Let <p be the normalized eigenfunction for g and put ip = a(p + <p' ( (<£> ><£>') =0), then 
<5l^ = 9o<fio yields 

oGl^o + Gup' = ag ip + StaV 5 '- 

Applying P and Qo, we have 

ag - a((p , D L (p ) - ((p Q , D L (p') = ag 
-aQ D L (p + Q G L (p' = g <p'. 
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Because of QqGlQo ^ QoGlQo ^ 9i < 9o, 9o — QoGlQo is positive, invertible and 

iff = -a[g - QoGlQo]~ 1 QqD l ipq, (3.11) 

go-go = (<Po, (D L - D L Q [g - QoG L QoY 1 QoD L )^o) 

= (Wi^il-WtQoigo-QoGLQoj-'QoW^Wtipo). (3.12) 

For brevity, we put 

X' = W* L Q [g - Q Q G L Q Q \- X Q Q W L , X = W* L Q [1 - Q Q G L Q Q \- X Q Q W L 



and 

Then we have 
and hence 



X = W* L Q [g - QoGlQ^QoWl. 
X -X' = -XX', 



X = X'(1 + X')- 1 and 1 - X = (1 + X')-\ (3.13) 
Together with W£ipo = Vl — e~fL~ d / 2 , we have 

go - g = L- d (Vl-e~f, (1 + X')~ Vl 1 e~f) (3.14) 

from (|3.12jl . Now, we want to replace X' by X in the right hand side. From 1 — go = 
0(L- d ), g -g x = 0(L~ 2 ), £^o V(l " gk? < K L ) OW, Lemma E^i)), it follows that 

\\X'-X\\ = (l-~go)\\W* L Qo[go-QoG L Qo]- 1 [l-QoG L Q }- 1 Q W L \\ 

= (l-~g ) sup (^WiQoigo-QoGLQoVil-QoGLQo^QoW^) 

ll*[[a=l 



< {!-%) sup ^|(<^Vl-e-/0)p 



9k 



k^O 



(9o - 9k)(l - 9k) 



^ (I -go) sup — > - 

\M\2=i Ld go - 9i f^o (1 ~ gk) 

= ||l-e- / ||iO(L- M ^(L)) = o(l). (3.15) 

Together with the similar estimate ||X|| ^ p c ||l — e~^||i(l + o(l)), we have \\X'\\ ^ p c ||l — 
e-^||i(l + o(l)). Thus (jHl yields 

' c-'lli 



L d (g - go) = (v 7 !^, (1 + X'y'Vl^) > ^-(1 + o(l)), 

l + p c \\l-e '||i 

which is the lower bound of (iii). The upper bound of (iii) is obvious. From 

iVl-e-/, (1 + XTWl-e-f) - (y/l-e-f, (1 + X)-Wl-e-f)\ 

^llv^^llllKi + ^iiiKi + xo^lliix-^i^oCi), 
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we get the first equality of (i). Replacing X' by X in (|3.14|) and tracing the argument back to 
(j3.12j) , we get the second one of (i) . (ii) is an immediate consequence of g ^ g and (j3.12|) . 
(iv) Clearly, a = (<p Q , (p ) = a'. As for (j3.13|) . we have 

(9o - QoG L QoY l QoW L = (g Q - Q Q G L Q y x Q Q W L {\ + X')' 1 . (3.16) 

This and estimates similar to ()3.15|) derive the bound 

||y?'|| 2 = a 2 ((p ,D L Q [g - QoG L Q ]~ 2 Q D L <p ) 

^ a 2 1 1 W* L <p Q 1 \l 1 1 WZQo [g - Q G L Q }- 2 Q W L | | 

= a 2 ||^o||2ll(l + ^ , )" 1 ^Qo^o-g G'Lgo]" 2 go^ L (l + X / )- 1 || 

= a 2 0(L~ 2 H(L)) 

from (|3.1H1 . Now the bound for 1 — a 2 is obvious. □ 
As in I, we use the generalized Vere- Jones' formula j^/J 1ST] in the form 



±- f per(J(x i ,x j ))» j = 1 \® N (dx 1 ---dx N ) = t _^Det(l - zJ)^ 

J J Sr (0) 



where r > satisfies \\rJ\\ < 1. S r (C) denotes the integration contour defined by the map 
9 1— > ( + r exp(i9), where 9 ranges from — tt to ir, r > and ( G C. Then we get 

eB r _ </ , e> - | _ z_l Det[l - z G L ] km D^l 1 - *G L {\ - gpg^ - l)]' 1 ^/2vr^+ 1 



Det[l - 5 G L ] / Si(Q) Det[l - z G L (l - ^G?^- 1 ^ - 1)] d v /2m V 



(3-17) 

The positive real numbers zq = zq(L,N) and zq = Zq(L,N) are chosen as the solutions of the 
equations 

Ti n [z G L (l - zoGl)- 1 ] = Ti n [z G L (l - ZoGl)' 1 ] = N. (3.18) 
In fact, the following lemma holds. 

Lemma 3.5 (i) z$ = Zq(L,N) G (0, 1) is uniquely determined by the equation 

Tt[z G l (1-z G l )- 1 ]=N. (3.19) 

(ii) zq = zq(L,N) £ (Q,% 1 (L) ) is uniquely determined by the equation 

Tr[z G L (l-z G L )- 1 ]=N. (3.20) 

(iii) ^ z - z = 0(L- d ) 

(iv) l-^ = (l + o(l))^(p-Pc)- 1 

Proof : Let H(zq) and H(zq) be the left-hand sides of fl3.19j) and ()3.2()j) , respectively. Since 
H is monotone increasing continuous function, H(0) = and if (1 — 0) = 00, (i) follows, (ii) 
is similar. The first inequality of (iii) is a consequence of H (z) ^ H(z). Before showing the 
second inequality, let us make the following remark on the thermodynamic limit ()3.5|) . 
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(a) If and only if p < p c , {z (L, N)} converges to z = z* G (0, 1), the unique solution of 

P = I J^y aiip;z) - 

(b) If and only if p > p c , L d {l — z ) — > l/(p — p c ), hence limz = 1. 

(c) If and only if p = p c , lim z = 1 and L d {l — z ) — > +oo. 

To show (a - c), note that 

z {L,N) f dp {L) 

= Tr[z (L,N)G L (l - z (L,N)G L )-']/L d = N/L d -, p. (3.21] 

We have that 

dp (L) x f dp 



(2n) d 1 ^ J Rd (2rcy 

for lim^o — z G [0, 1] by the dominated convergence theorem, and that the limit is a strictly 
increasing function of z. ( See Lemma f3. 21 ) If lim^ — z * £ [0, 1), the limit of (|3.21j) tends to 

If lim^o = 1, then p = p c + \imzo/L d (l — zq) ^ p c . Now suppose {z (L, N)} does not converge. 
Then by taking converging subsequences having different limits, we deduce a contradiction to 
(13.21)1 . Thus we get the classification (a - c) and (iv). 

Now we have the second part of (iii) using Lemma l3~4T iii) . 

2 = 1 - 0{L- d ) ^z < ~g Q l = 1 + 0(L- d ). □ 



In order to understand the subsequent arguments, it is helpful to keep the followings in 
mind: 

A, = 1, 9i = l- 0(L- 2 ) > QoG L Qo > QoG L Qo (see §J^) 
g = 1 - 0(L- d ) (Lemma E3t iii)) 
Zq = 1- 6{L- d ) z = z + 0(L- d ) fLemmaECTiii. iv)) 
bf\D LV V)=g?>\\l-e-%fL* 



Lemma 3.6 



(i) P [l - ZoGl^Po = (-^r-r- + O(L- d £(L)))p , 

(ii) ||Q [1 - z G L ]- l \\ = \\[1- z G L ]- l Qo\\ = 0{^I{L)), 
\\Q [1 - z G L ]- l \\ = ||[1 - z G L ]- l Q \\ = O(V^I)), 

(iii) Tr(Q [l - z Q G L \- l D L [l - z G L ]~ l Q ) = 0(L- d £(L)). 
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Proof: (i) By lemma HOliv) , we have 
\(<p , (1 - 5 Gl)~Vo) - (1 - zogoT 1 ] = \(a<Po + 0', (1 - z Q G l)' 1 {aip G + iff)) - (1 - z Q g Q y l \ 

^ ^—S^ + W, (1 - *>G L )-*<ff)\ ^ ((1 - ^o)- 1 + (1 - ~zmr l )0{L-*%L)) = 0(L~ d i(L)), 
where we have used 

X — + : \ - < 2 + Tr[f G L (l - Sq^)" 1 ] = 2 + N = 0(L d ), 



1 - z g 1 - £ <?i 
in the last step. 

(ii) Note that Qofio = V 9 ' i n the notation of lemma E31^iv). Then we get 

1 



||Q (1 - zoGlY^ ^ + — = 0(L d y/L-^e(L)) + 0(L 2 

1 - 2 O 0o 1 - z gi 

The second bound is obtained similarly. 

(iii) From the above remark, the left-hand side equals 



-E 



which yields the righthand side by Lemma f3. 3f i) . □ 
We need a finer estimate than Lemma f3.5f iii). 

Lemma 3.7 



(i) z -z Q = (l-g )(l + o{l)), 

(ii) 1 - z g' = (1 - z g )(l + o(l)) = (1 - zq)(1 + o(l)) 



l + o(l) 



where g' = 1 -((p^\ D L <p^) + z ((pf\ D L Q [1 - zqQqGlQoI^QqDup^). 



Proof: (i) Let us begin with 

= N - N = Tr[5 G L (l - ~z a G L )~ l - z G L (l - z G L y 1 } 

= (tpo, ((1 - zoGl)- 1 - (1 - zoGl)- 1 )^) + Tr[Q ((l - 5 (5 L ) _1 - (1 - z Q G L y l )Q Q ] 
+Tr[Q ((l - zoGl)- 1 - (1 - ZoGy^Qo]- 
The first term of the right hand side equals 

(1 - ~z ~g o r ~ (1 - Wo)' 1 + G{L-H{L)) = %^^~ g 'f + 0{L~*l{L)) 

by Lemma f3. . On the other hand, the second term has the bound 

(So - ^o)|Tr[g (l - z Q G L )- l G L (l - z G L )- l Q ]\ 
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< Z -^^\\~ Zq G l {1 - zoG L )^|| T ||(l - zoG^QoW = 0(L~ d L d ^I(L)) = o(L d ) 

by Lemma [3.5f iii. ii) and Lemma [3.6f ii). The third term can be estimated as 

|Tr[Q ((l - zoGl)- 1 - (1 - z G L y l )Q }\ = ^ |Tr[Q (l - z G L )- l W L Wi{l - z G L )- l Q }\ 

= ^o|Tr[Q (l - z G L )- l W L {l + z W[(l - zqGl^Wl^WKI - z G L )- l Qo]\ 
*C *b||Q„(l - z G L )~ l W L Wl{\ - z G L )- l Q,]\\ T = 0(L- d £(L)) = o(L d ), 
where we have used a equality similar to (j3.16|) and Lemma fe^f iii). Thus we have 

Zpjgo — go) — (zp — Z )go _ 

(1 - z g )(l - z g ) 

On the other hand, (1 — zogo)(l — z g ) = 0(L~ 2d ) holds. Thus we have 

za{go - go) - (zq - z )g = o(L~ d ). 

Note that go — go is exactly of order L~ d by Lemma E3fiii), we get the desired estimate, 
(ii) From (j3.12j) . we have 

\go-g' \ = K^DLQoKgo-QoGLQor'-Cz^-QoGLQor^QoD^o)] 

= \(f , D L Q (g ~ QoG L Qo)- 1/2 [Cz 1 - g ){% 1 - QoGlQoY 1 ]^ - Q GlQo)- 1/2 QoD lVo )\ 

< l^- 1 - g \ 1 1 (z 1 - QqGlQoY 1 1 1 fao, D L Q (g - QoGlQ^QoD^o) 

^ 0{L- d )0{L 2 ){^,D L ^) = 0(L 2 - 2d ) = o(L~ d ), 

where Lemma I3.4f n) has been used in the last inequality. Hence, we obtain 1 — zog' = 1 — 
Zogo + o(L~ d ). On the other hand, we have 

1 - zogo = l — zo+ [2q(1 - g ) - {zo - z )} 
L d {p - pc) 

thanks to Lemma f3.5f iv) and (i) above. □ 
Put 

(AO _ z (L,N)gj(L) „ (A r) _ z (L,N)gj(L) 

Pi 1 . I T HT\ . I T \ ' Pi 



1 l-z (L,N) 9j (Ly ^ l-z (L,N)~ 9j (Ly 
then we have Yl'jLoPf*^ = Yl'jLoP^ = N by Lemma l3~HT i. ii), 

pi N) =0(L d ), p^ = d(L% ^7^ = 1 + 0(1) (3.22) 
by Lemma f3.7f ii) and 

A N) =0{L 2 )>A N) >---, P{ N) = 0(L>) >p^>.--. 
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Lemma 3.8 

1 dr] l + o(l) 

, l(0) Det[l - z (L,N)G L (l - z (L,N)G L )^(r] - 1)] 2 7 rir 1 "+ 1 ep f) 

1 dr] l + o(l) 



5l(0 ) Det[l - z (L,N)G L (l - z (L,N)G L )-i( V - 1)] 2m V ^ e pf) 

Proof : Set R.W = R.W = L^ 2 )/ 2 . Since E^pf^l + pf°) = Tr[z Q G L Q (l 
zqQqGlQq)- 2 } < EJLi W( x - #j) 2 > w e get 

pf)(l+pf)) 



Po 



by = 0(L d ) and Lemma I3.3IT J. Then Lemma A. 2 yields 

1 cZ?7 l + o(l) 



the l.h.s. of the 1st eq. 



For the second equality, we notice that pj ^ (1 + o(l))pJ- holds for all j = 1,2, ••• 



because of z , z = 1 + 0(L~ d ) and < gf ] ^ 1 - 0(L~ 2 ). Together with lET22|) . we have 



'.7 



^ E - =1 pf) ( i+pf)) ^^^Er^pf^+pf^ 



< + 



~(7V)2 ^ V 1 ^ " (AT) 2 

Po Po 
Thus the second equality also follows from Lemma A. 2. □ 

Now we have 

Det[l - ^ G L J 

from (|H.17|) . ()3.22|) and the above lemma. Since Po,Qo and Gl commute, Det[l — z Gl] = 
(1 — zo)Det[l — zqQoGlQo]- We use the Feshbach formula to get 

Det[l - z G L ] = Det ( P ° ~ z ° P 9, GlP ° ~ z o p o g lQo 
— ^^QqHl [Qo — ZqQoGlQo] 



xBetp n L [P - z P G L P - z PoG L Q {Qo - z Q G L Q ) z QoG L P { 



o 



Det[l - ZqQoGlQo] 



x (1 - z [l - & ( L \ D LV ( L) ) + z (v ( o L \ D L Qo[l - zoQoGlQo^QoDl^)]) 

where Det is the Fredholm determinant for operators on Tii, and Detg -^ L for operators on the 
subspace QqHl etc. Now from Lemma l3~7T ii) and Lemma EToT iii. iv), we get 

EfJe"<«>l = Z -j P--*)^ 1 - ZqQoGlQo] (1 + o(1)) 
' L J ~4 (1 - % )Det[l - zoQoG L Qo] 
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ztf Det 1 - zqQoG l Qo \ f m * 

■^Nr—— — — ~ (1 + o 1 

4 Det 1 - zoQqGlQo] 



( ^ ~ z ° at , m ^ Det t 1 - ^oQoGlQo] Det[l - Q G L Qo] Det[l - Q G L Q ] 
exp iv+o 1 — — F — — — — : = = . (3.23) 

V zo " Det[l - g G L Q ] Det[l - Q G L Q ]Det[l - z Q G L Q ] 



Lemma 3.9 



m Det[l - zqQqGlQq} _ ( l^zo , N (aQn , m \ 

W Det[l - g GT L Q ] " eXp ^ z {N Po ) + 0{1) > 

Det[l - Q G L Q J 

Det[l -Q GxQo] n , M , „ U1 , 
(m) _ [ _^-^ I = Det[l + ^](l + o(l)) 



Proof : Put = — log Det (1 — zQ GlQ ) = — X]^=i l°g(l — z<7j), an d we have 
where zo £ (20, !)• Hence we get (i) by 

9j N- p 



and 



*"<W - *)' - 1 f^S < I fr^ = owl)) = o(i), 

where Lemma lH.Hf i) has been used. Similar argument and ^ Sjli 

Pi) y ield (h)- 

(iii) Thanks to the product and cyclic properties of the Fredholm determinant, we have 



Det[l - QqGlQo] 



T>et[l+Q Q (G l -G l )Qq{\-QqG lQqY 1 } = T>et[l+WlQ Q {l-Q G L Q )-*Q Q W L ] 



Det[l-QoG L Qo] 

= Det[l + Vl-e-fQ G L Q (l - QoGlQ^ 1 y/l - e~f]. 

Note that L 2 (A L ) can be identified with an closed subspace of L 2 (M. d ) naturally. By this 
identification, we regard Gl and y/l — as operators on L 2 (M. d ). Now for (iii), it is enough 
to prove 

A L = Vl-e-fQ G L Q (l - Q G L Q )-Wl-e-f — > K f 

in the trace norm. In the following, we show Al — > Kf strongly and ||Al||t — > \ \Kf\\T- Then 
the Griim's convergence theorem [Slj yields the above. 
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For ip,(f) e L 2 (M, d ), we have 



\W,(A L -K f )<f>)\ 



dx / dyip(x)y/l- e~f^(j)(y)Vl - erM 



([QoG L Q (l - QoGlQo)- 1 ]^^) - K(x,y)) 



^(|^|,V / l-e-/)( V / l-e-/, 



sup \[QoG L Q (l-Q G L Q )- 1 ](x,y)-K(x,y)\ (3.24) 

x,j/£supp/ 



< ll^lbl^lbllVl-e^ll^ sup \[Q G L Q (l-Q G L Q )- 1 ](x,y)-K(x,y)\, 

x,y£suppf 

which tends to 0, by Lemma l3.3n v). Thus the strong (in fact the norm) convergence has 
been proved. For the convergence of the trace norm, we use Lemma I3.3f iv) again and positive 
self-adjointness of operators A L and Kj to get 

\\A L \\ T -\\K f \\ T = Tr[A L -K f ) 
= [ dx(l-e- f{x) )([Q G L Q (l-QoG L Qo)- 1 ](x,x)-K(x,x))^0. 

□ 



By Lemma [3.51 (iii, iv), we have 

— Zq 1 — Zq 1 — 5q [Zq — Zo)(l — Zq) 



-N + 



-N 



-N 



1 



Zq Zq Zq 

Applying Lemma f3.4f i) and Lemma r3.7f ii) to the righthand side of 



N = o(—). 

z z y N J 



I- Zq 1 - Zo „ I -90 

-Pa H ~ — Po 



Zq 



zq 



1 - z g 



we get the formula 



K N [e- <M> ] 



exp ( - (p - p c )(Vl - e-f, [1 + W[Q (l - QoGlQ^QoWl}- 1 ^!^^) + o(l)) 

Det[l + K/] 



(3.25) 



From the convergence WIQ (1 — QqGlQq) 1 Q$Wl = A L — > Kf in the thermodynamic limit, 
we have proved the theorem. 
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A Complex integrals 

Lemma A.l For ^ x ^ 1 and p ^ satisfying ^ px < 1, we have 



^( 1+ ^ 1 -p,)>exp(- ^ 2 ^y) ^). 

Proof: Put /(a;) = log(l + xf(l - px), then 

f'M = _ g f'M = P(l+P)(l+P^ 2 ) 

71 J 1 + s 1-prr' /W (l+x) 2 (l-px) 2 

hold. So we have /(0) = 0, /'(0) = and > /"(0a;) ^ -p(l + p)(l + p:r 2 )/(l - px) 2 for 
G (0, 1), which imply the result. □ 

Lemma A. 2 Let the collection of numbers {pj N ^ }j,N satisfies 

oo 

P ( o N) > p[ N) > pf ) > ■ ■ ■> Pf ] > ■ ■ ■ > 0, £ pf> = AT. 

i=o 

Suppose that there exist a sequence { R^}n^ and c E (0, 1) snc/i that 

1 < i2<"> < cp^ (l A 4)) ' Jin^V^e^' = 



Pi 



and 



T/ien 



holds. 



y ^E^^d+Pf) n . , I0g(l+C) 

hm — = 0, where c = - 

N^oo Po C 



lim (h 

J s 



dr) 



5l(0) 2«^ + i n oo o(1 _ p w (r? _ 1)) e 



Proof: We omit the superscript (TV) here. 

Note that po — > oo and i? — > oo as N — > oo. 
By the preceding lemma, 



1 > A [(1 + £)" (1 - %)] .e XP (-|Ml±M (l + ^)J (l _ ^ n 



J'=l 

So the assumption on i? implies 



n[(i + «r(i_^) 

jLi LV Po' v Po ' 



AT->oo 
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Similarly, we have 



£=i L Po' v Po 



1. 



Now let us deform the integration contour of the complex variable r\ to two parts 

= - I + I = h + h- 

Si(o) Js (R _ 1)/PO (i+i/po) Js 1+R/PO (o) 

Ii is obtained by the residue at r\ = 1 + l/po : 



h = -Po 



Po 



1 + - (-Pb)lK 1 - 



Po'J 



= 1 



1 \ -po-l 



Po 



J-J- LV p / V p 



e" 1 . 



Jo can be estimated as 



\h\ < 



-£In[(-£)1^((-£) 



le" - 1 



sM'+-r"ii-«r'[n( 



i + 



Po' L " V P0 / V PO 

since (1 + R/p ) Po > (1 + c) R//c = e c ' fi and the assumption. 



AT— >oo 



o, 



□ 
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